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We report on calculations of the hadronic coupling constants (jp™- and qnntt based on lattice QCD with four 
flavors of dynamical staggered fermions. By computing 2-point and 3-point Green's functions we have been able 
to determine these coupling constants ab initio from QCD; the results gp-^ = 4.2 ± 1.8 and qnntt = 14.8 ± 6 are 
compatible with the experimental values. 
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Recent lattice QCD calculations of the had- 
ronic mass spectrum have given encouraging re- 
sults, reproducing many of the experimental val- 
ues even at a quantitative level. This motivates 
the launching of more complicated analyses such 
as the calculation of properties of the phenomeno- 
logical hadronic interactions. The aim of our 
work is a calculation of coupling constants which 
parameterize the strength of the interaction be- 
tween hadrons at low energies. 

We have based the analysis on the staggered 
fermion action with four degenerate flavors of dy- 
namical fermions. The lattice size is 16"^ x 24 
and we have generated so far 85 configurations at 
/3 = 5.35, m = 0.01 using the hybrid Monte Carlo 
algorithm. Details with respect to the production 
of configurations can be recalled from Rcf. [Q . 

Hadronic Couplings in the Continuum 

Interactions between hadrons at low energies 
can be described by effective Lagrangians in 
which the hadrons enter as the "elementary par- 
ticles" of the theory. The strength of the inter- 
action is parameterized by effective coupling con- 
stants whose values have been determined by ex- 
periment. In our case the effective interaction 
Lagrangians are 
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The Lagrangians are scalars with respect to 
Lorentz and flavor transformations, the latter be- 



ing achieved by contracting the flavor indices with 
the appropriate SU{Nf) coefficients Cabc and c'^i,^. 

The (7pmr coupling can be fixed by measuring 
the width of the p decaying into two pions giving 
the experimental value g^^ ~ 6.08. 

The coupling constant gNN-n can be extracted 
from a low energy phase shift analysis of ttN elas- 
tic scattering, leading to the value of 5^^ ~ 13.4. 

Hadronic Couplings and 3-point Functions 

N-point Greens function are calculated on the 
lattice by means of the path integral method; the 
hadronic coupling constants demand the calcula- 
tion of 3-point functions. 

To give an example, we consider the case of the 
pTTTT coupling. Here the following 3-point function 
must be computed: 
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e^ is the polarization vector of the p meson and 
the tilde indicates that the corresponding field op- 
erator is defined in momentum space. Inserting 
complete sets of particle states, using the LSZ for- 
malism and applying the phenomenological inter- 
action Lagrangian, the following relation between 
Fs and gpT„^ can be derived (ti '3> t2 ^0): 
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Two terms in this formula depend on the specific 
process: The tensor Cabc and the kinematical fac- 
tor 2\qp\^. The other terms are analogous for 
the calculation of the NNtt coupling constant. 

To determine the value of gpmr the knowledge of 
the masses rrii, the energies Ei and the renormal- 
ization constants Zi of the particles is required in 
addition. These observables have been evaluated 
in a previous analysis of the spectrum by com- 
puting the proper 2-point functions ||l| . 

Hadronic Operators on the Lattice 

Different approaches can be applied to the 
construction of the appropriate hadron opera- 
tors from the single-component staggered fermion 
fields X- The application of group theoretical 
methods is an adequate and mathematically rigid 
method to construct hadron operators with stag- 
gered fermion fields. In our case, however, a more 
intuitive way is preferable to simplify the calcula- 
tion of the group coefficients. First, quark spinor 
operators are constructed which are then used to 
build up the hadron operators. In this way the 
flavor of the lattice operators can be defined in 
terms of the usual isospin, strangeness and charm 
quantum numbers, allowing us to use standard 
tables of S'[/(4) Clebsch-Gordan coefficients ||J|] 
to calculate the group tensors. 

The quark fields are defined on a coarse lattice 
whose points y consist of hypercubes of the orig- 
inal lattice and the "/ine degrees of freedom" are 
used to provide them with spin and flavor indices: 
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Here T^" 



(7i''''74*) with a, a denoting 
spin and flavor indices of the quark, respectively. 
Mesons are of the form 
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where the multi-indexed products of 7 matrices 
Ts and Tp determine the spin and flavor content 
of the state. In terms of the staggered fermion 
fields X these operators are non-local if the spin 
multi-index S and the flavor multi-index F are 
not equal. 

The local baryon operator which corresponds 
to the nucleon in the continuum limit can be con- 



structed as follows 

B^^iy) = ^e,,kJ2(TAq^{y)r\ 

A 
x{qj{y)[CTA®C*n]qk{y)}. (6) 

C denotes the charge conjugation operator and a 
summation over color indices i, j, k is understood. 
a and a are spin and flavor indices of the baryon. 

Hadron Masses 

We merely recall here the results of the analysis 
for the spectrum of the required particles; for de- 
tails see y . A few remarks, however, ought to be 
added: (i) Some of the results are technically new 
as the masses and energies had to be computed for 
polarized p mesons and nucleons. (ii) The parti- 
cle energies have been determined for states with 
minimal lattice momentumR qmin — ^ ^^ the z 
direction. Fixing the average p mass to 770 GeV, 
this corresponds to qmin ~ 550 GeV in physical 
units, (iii) For the p we used a local-time (LT) 
as well as a non-local-time (NLT) operator, the 
latter supressing the unwanted parity partner of 
the p. In the case of the nucleon we used a wall 
operator for the source. 





m 


E{q) 


n[2] 
n[7] 


0.271(5) 
0.346(15) 


0.50(3) 
0.53(3) 




0.55(5) 
0.53(2) 


0.70(8) 
0.65(9) 


Npoi 


0.76(7) 


0.84(8) 



Table 1: Masses and energies in lattice units. The 
numbers in square brackets refer to Ref.||l[. 

The gpTTK Coupling Constant 

We have calculated the 3-point function eq. (g) 
with the TT meson being fixed at the time slice 
tjr = 5 for all values of tp. The p meson and one 
TT meson were given a momentum qmin in the z 
direction. To avoid statistical noise, the operators 
were chosen as local as possible. We started with 
a spatially local p meson and the local "lattice 
Goldstone"-pion, the latter having the additional 
advantage of being the lightest pion and having 



16 is the spatial extent of the lattice. 



no parity partner. However, due to the flavor 
symmetry the p cannot decay into two local tt's 
so that the second pion had to be chosen non- 
local. The calculation was performed by using 
both LT and NLT p operators. 

The results which are achieved by doing x" fits 
of the data to eq. (||) are: 



g^^ - 3.6 ±2.5 (LT) 
gpr^ = 4.7 ±2.8 (NLT) 
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Though the values are systematically low, within 
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Figure 1. The data and a x^ fit for gpT^{LT) 

the rather large errors they are nevertheless com- 
patible with the experimental value. Moreover, 
one should bear in mind that only one pion was 
light in the p decay so that the particles had to 
be off-shell. 

The gNN-n Coupling Constant 

We have computed the 3-point function using 
a polarized nucleon, a local lattice GS-tt and a 
wall-source nucleon operator (only half of the con- 
figurations have been analyzed so far). The wall 
source projects onto a nucleon at rest while the 
sink nucleon and pion have been given a momen- 
tum qmin in the z direction. The sink nucleon has 
been fixed at the timeslice t^ = S> and we have 
evaluated the correlation function for all times t^r 
of the pion. 

Unlike the pmr coupling, gNN^^ could not be de- 
termined yet directly from the data evaluated so 
far. This is due to the fact that there are two 



different couplings gi and g^ for the N — > Ntt 
process in SU{A). The result obtained from this 
analysis, is therefore a linear superposition of two 
couplings and additional information is needed to 
disentangle them. This information can be ex- 
tracted from the lattice^; however, for the time 
being we have used the experimentally known 
value of Fa/ Da = 0.59(5) Q to process the data. 
The final result for the NNtt coupling is given 
by: 

gNm = 14.8 ± 6 

The central value is much closer to the experi- 
mental data than the result for gpmi-; this may be 
due to the fact that the gNN-rr analysis does not in- 
volve "non-Goldstone" pions. However, the error 
is still large and must be lowered by increasing the 
statistics. Nevertheless, this ab initio calculation 
of the NNtt coupling from QCD is an exciting 
first step — we think. 
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Figure 2. The data and a x^ fit for gNN-rr 
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Analysis in progress. 



